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SUMMARY 


The present report is concerned with a theoretical study of the 
responee of a linear tw terminal network to a frequency puise modulated 
wave, if this network is followed by an ideal discriminator. We have assumed 
the modulating function to be of the forn 


where t is the time variable, and » a parameter characterizing the "duration 
of the pulse®. 


In Section I the response r(t) of a general two terminal retwork 
ia found for the special case when A 0/2 is an integer, where A @, ie 
the deviation from the carrier frequency @,° 


In Section II the output of the discriminator r,(t) is determined, 
under the assumption that A w/a: = 1. It ie shown that such hypothesis 
corresponds to cases of practical interest. 


Section III is concerned with oxpansions of the exponential inte- 
gral of complex argumert. This fimction enters in the expressions of r(t) 


In Section IV tho formulae obtained in Section II are specialized 
for the case of a single parallel RLC circuit taned on the carrier frequency 
mo, . It i¢ assumed that 2m >> and 2 w, >>vw, whare aw is the total band- 

dth of the circuls between half power points, 


In Seotion V asymptotic expressions are de.ived for the cases of 
wip —>roo and w/p— 0. For w/u—Po it is found that vhat time dependence 
of v,(t) approaches 8(t) as, of course, it should be. In such a case, the 
prrallel tuned circuit degenerates in to a pure resistance (w= 0), if it is 
assumed p 4 0. In deriving these asymptotic expressions the asseumpticns pre- 
viously made are of course maintained. For instance, it is always assumed 
a0, >>. 


Section VI is concerned with numerical computations, based on the 
fermulae obtained in Section IV. Oalculations have been carried out for 
-4&< ut <4, and for w/p—>o, 10, 4, 2, 1, 0.4, 0.1, 0.01 —>0. These cal- 
culations show that the time dependence of v,\t) is very similar to that of 
the modulating function 8(t), provided w/u >1l. If w/w <1 the input signal 
becomes very distorted. We can say that for o >wie > 1 the output 1s similar 
to the input, 1.0. 1/l+(ut)®. For 0.1 < w/y <1 the output is very critically 
&spenu2nt on time, and there:ore difficult to calculate. In the interval 
0 <w/p< 0.1 the output does not vary too much and it is similar to that given 
by the esymptotic formula for w/p-—90, 1.9. wt/ts(pt)e. 
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Introduction 


Transients solutions for linear networks, subject to a frequency 
modulated driving function, have been obtained for a few typical cases. 
Salinge considere the reeponse of an ideal band page filter followed by 
an ideal discriminator, to a frequency "step". Gold” extgnded Salinger's 
solution for the case of a frequency "squzre pulse®. Hok and Ekstein and 
Schiffman’ have investigated the response of a linear network, without dis- 
criminator, to a linearly increasing frequency. 


The solutions recalled were possible in virtue of the simplicity 
of the modulating function chosen. Salinger picked the sinplest type of 
modulating function, i.e. a constant (in Heaviside’s sense). Eok considered the 
next function in order of complexity, i.e. the Linear function. 


It seems therefore natural to ask: Is there gome other type of 
simple modulating function which, besides being of practical interest, will 
lead to integrations which can be handled? If one considers the formal ex- 
pression for the response of a general linear network subject to a frequency 
modulated driving function, it is easy to see that modulating functions of 
the form 
S(t) = 1 i 


t+ pe te : 


where t represents time and pp a parameter, lead to exnressions which can be 
integrated, at least for values of the frequency swing integral multiples 
of &. The function S(t) is aleo suitable to represent a “pulse”, which is a t 
type of modulation of practical importance. We have therefore chosen this 
function for our investigatior. 


i H. Salinger, "Transients in Frequency Modulation", Proc. I.R.B., August, 


19a. 


B. Gold, "Transients in Frequency Modulation", Report R-154-47, PIB-103, 
Microwave Research Inetitute cf the Polytechnic Institute of 
Brooklyn, 1947. 

3 G. Hok, "Response of Linaar Resonant Systems to Excitation of a Frequency 

Varying Linearly with Time", Jour. App. Phys., March 1948. 


H. Bkstein and T. Schiffman, "Responses of a Linear Network with a Linearly 
Variable Frequency as obtained in Sweep Frequency Test", Pro- 
ceedings of the National Electronics Conference, Yol. VII, 1951. 
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Section I - Response of a Two Terminal Network 


The response of a two terminal network can be expressed, by means 
of the superposition integral, as follows 


t 
r(t) = i fie) gt ~ T) at (1) 
"00 


vhera: 
r(t+) is response 


f(t) ie driving function 


e(t) is the Green's function, 1.e. the response of the netwrk in 
question to a unit impulse. The lower limit of the integral has been assumed 
to be ~ 00 since the function f(t) that we shall consider starts at t= ~@. 


For a non-impuleive network the Laplace transform of g(t), that we 
shall indicate with G(p), is a proper rational fraction, and can be expanded 
as follows: 


n A, 
te) ey peo (2) 
sal 


where p_ are the roots of the denominator of G(p), which are assumed to de all 
distinct. We now assume a complex driving function of the form 


a(t (wt+A 8(t) dat) 
f(t) «= x Ri = F = “ot 5 we (3) 


The instantaneous frequency, ies defined uniquely, since es is constant with 
time, as 


° 


ag 


at — ®, + A a S(t) (4) 
If now we assume 
S(t) = . (5) 
1+ (ut)® 


the frequency of the driving function will vasy {= tne pulsc like m*nner de- 
picted in Fig. 1. 
5(t) 
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t 


Since We 8{t) at = : (tan7? 


pt + =) Eq. (3) can de written: 


Me t+ SS (tan7? pt +2) J 
f(t) = Foe bie . e (6) 


If now the function (3) is apvlied to the two terminals of the 
eysten (2), the typical » comporent of the response can be written, by usine 
the superposition integral, as follows 


= P,* ie a aa 
Be A, e ay @ t(c) ae 
and using Eq. (6) Ao 
Pe pt pt +[ t+ Jj R 2 ten? pe 
yr, = Te A. c) = 6 av (7) 
Where we have let 
Aw 
3 a 
— ae fo tee Fo Te (8) 
8 0 8 2 n) Mra 
4f now we recall that 
-1 aie 1+j2 
tan 2 xz 2g log Ley e 
Eq. (7) becomes Ate 
o ?,¢ s Pee are = 2 
ro ay Ale aes (pts apr ot) d (9) 


The integral appearing in Ec. (9) can be easily evaluated 1f Aw /a is an 
integer. If we let A wo, /2u = @, we are concerned with the following integral 


qa fit 
le Get) 0 * at (10) 


To tranaform this integral we notice that, in general, we can write 


(z+ 1)%7 a E, Hi 
= + a Hane a (12) 
(zg -1)% Z-1 (Z - 1) (% - 1) 
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where 
q-r = 
a ea rel,2,...4 (12) 
(q-r)t a2 tl 


since the function on the left side of Bq. (11) has a pole of order q at Z=1l. 
Yorforming the differentiations indicated in Eq. (12) we get 
q-1, ,r-1 
‘= CG - : 7 


¥q.(10) therefore becomes 


i (Leis) A dt = (-1)% i (Luts 2)* sist 


-o l- j ut -o jut-1 ae 


t [e-a te seu 
- (a)? f oN, aH Lt EH. ... 
i ted ea (t+ dé 3 (3 n)* (ds t)? 


t- a r+ 
a e © at 
(J w) (4 ¢)t 
Letting now 
zetet, 2- Bo s- dt bate 4d, ataas 
we get 
P on, + 8 
‘ ya 7 r ae aa” pt 
Si, GEE eH at 2 (rte + oo NS = ad ju Zz 
w+ 8 2 an i 
ah eA tet “El 5) a * a 
z (3) 3 (Ju) & 
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where the path of integratioxr in the complex Z plan is indicated in Pig. 2. 
j vs 


Z 2 LAME 
4h 


e 
7 


REAL Axis 


t 


Pig. 2 


All the integrals appearing on tha right hand side of Eq. (13) can 
de expressed in terms of the exponential integral 


(j/o+t) [8 


e 
- i tp i Ae c ae —zj- «Oa (14) 


this is apparent if we recall that in general 
2 
os 


gz 


z 


ri 4 1 
Na = 8 mi gm-l * (m-1)(m-2) gm-2* °° * (m-i)! 


r m-l [2 
a Geel) 77 a m>l- (15) 


Bq. (15) 48 obtained by successive integration by parts. In applying this 
formula to Bq. (13) we have to recall that 


Om eee fe a, [-», | >° 


and therefore 


e * — 0 when B [z]—»-o . 


; 
| 


oS 
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By inserting (15) into (13) we finally attain the following 
expression for (9) 


ave 
Bide ark Pa a rar ho eet i case BSS | 
8 2 8 yt (iiayog == 
q-1 
(ju)? (q-1)! . 
Cee Ta ey 
Pr, (sw)? tad (5m)? % (eo dy? 
Pe Rose Rae Ree! Oe COIN VELL ae nee aa 
f (t+2 (ju)t od (t+ dja? (q-1)(a-2) (t+ ayer? 
r a2 
Ske ee ak ) (16) 
(q-1)$ tea 


Equation (16) shows that the response of a general network to a 
frequency modulated driving function, whose modulating function is of the form 
(5), can be expressed in terms of exponential integrals of complex argument, 
provided the ratio A w/e is a positive integer. 


a 
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Section II - Output of the Limiter-Discriminator 


Equation (16) gives the sth component of the response. What is 
usually required, however, is the output of an ideal limiter-discriminator 
network, to which the response of the system is applied. For ideal limiter 
discriminator we mean a device, which first reduces the amplitude of the 
applied signal to a value constant with time (limiter), and whose output is 
proportional to the time derivative of the phase of this constant amplitude 
signal (discriminator). 


Formally it is a simple matter to find such derivative, but to carry 
out the actual calculations is usually very cumbersone. 


In general, if we write the response as 


He = J Re So ese) = pest (17) 
‘=1 e=1 
its phase angle can be written 
+ = tn? 3 (18) 
and therefore 
rit. x oo = a(t) (19) 


where K is a constant which is dependent upon the characteristics of the dis- 
orininator. 


From the preceding section it is apparent that to derive a general 
expression for the responses rg for any integer q is e very complicated matter. 
In what follows we shall therefore only consider cases where 


A®, 
2h 


That such an asswumption refers to cases which can be of a practical 
importance can be easily shown. 


q = a2 (20) 


From Eq. (5) it is evident that we have to define the duration of 
the pulse. A reasonable definition for this quantity can be the following: 
The duration of the pulse is defined as twice the time T necessary to reduce 
8(t) from 1 to 1/10. Using thig definition we can relate » with the duration 
of the pulse 2T. We have 


0.1 = 2 
1+ (ur)? 


B = (21) 


or 


| 
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Yor a pulse duration of 1 » sec (2T = 10° sec), we get from Iq. (21) 

b= 6.10° 

and taking into account Eq. (20) 
A», = 12.10° rad/sec. 


In other words with a pulse duration of 1 » sec, to fulfill condition (20) 
the frequency deviation must be about 1.9 Mc/sec. 


Asewmption (20) simplifies very much expression (16), end ve have 


F,t- d icc Se ls G ov 
Foe -a te bes 2 ] (22) 
a 
since according to the second of (8) Fo = P, eo” az t,: If now we let 
J 
Bh + iP. ; A= |ale ° (23) 
and recall the first of (8) and (14), we get 
Mo Ms <s 
= , —,° + + j g t+ 6 + wv 
a T, la, | a(t. + J The) d - ‘ ° 7 
. ~<,+3 (o, - @,) 4 +o, | (4) 


~, + (w, ~ py 


The second term of this equation is easily recognized to be the steady state. 
Separating the real and imaginary part of the right hand side of Bq. (24) we 


% 
- t at - b, 


ORiae car eee 


(25) 


R-HI“-H2, PIB-238 9 


o -B 
moh es ° 7 8 
xX = F, la, a [ - Ie com Y + a ein vi 


ee 1 rire [(e, - f,) cos 5, + o<, sin 5,| 
oy * (w, f,) 
(26) 
where 


Yo B,t+e+ = and J d+, + (27) 


the derivatives of RB. and x are as follows: 


m4 


=<, 
a = § ! 
Ri r, lale At st Yetti 008 Y +41, 008 ee sin | 


~,* r |a 
+0 [Fy 08 v tT, sein Y,| + ites [~pegein 5,0, (2,-A)o085, | 
so fs 


(28) 
a,-P, 


e a | 
As] elt t 
x = t, lala c) yt, gin v; oe Yeth ty goes Yrt sin LP 


Fr fa| 
+ =< [-t,,.28 voy ge | +t (apf, ein =~», 2084 | 


~t, He,- f,)° 
(29) 


Formulae (25), (26), (28), (29), together with (19), allow the calculation of the 
output of the limiter-discriminator, when the network in question is driven 
by the function (€), and Ao,/as 1. 


EE TED TT ae Fey a nST EDD SERNA PTE SRR Fy: 1 TE RRR TaenP ENED Wnnerterrere steerer. IR AERENNPTRARG SwT AL RENO  NT APPL 
e 
o 


eee. ee een eee ES a PS SA ees EE Et wr eee ee 
= — _ 


ve 
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Section III - Expansions for the Inte 14 
The integral (14) oan be easily transformed into a standard form as 
followe: 
S/wot _t2 hy 
-n+ (+00 a) ic 
where » 
-P (det) = fe slo, - 27] Coe 8) (32) 


oN + j¥ 2 Ge 


wo - 
Meat + é Hs =-(, -@) + (32) 
Yow we have 


ry - () 2 
T s C+ log a-f ee = O+logA- Cu — (33) 


where C = 0.5772 is the Buler-Mascheroni constant. 
Since the path of integration does not encircle the origin (Fig. 2) 
logh = leg +Jée (34) 


we can therefore write 


Ie I, + jt, = 0+ log (+ uy eeu Satis 


(35) 


A useful asymptotic expansion for the integral (30) ie easily obtained by suc- 
cessive integration by parts, and it is the following 


-d 
e oS ie 31 
Aa ae Sonate - 3 Paes (36) 


FF PR OPP EE FO ER FTE DT TF OT EET 
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It is easy to derive an expression for I' which appears in Bq. (28) 
and Eq. (29) 


=A - J (w, - B) ace 
rn,42,4! $b &— [-<-s(u,-6)] = p 


at aa «at 


+J 


ee 
ea] 


(37) 


Section IV - Single Tuned Circuit 


We want now to apply the formula of the preceding section to cal- 
culate the output of the discriminator, for the case sketched in Fig. 3. 


(th Yay) 


LUNMTER OLSCRMIMINATOR 


Fig. 3 
We shall assime that the circuit is tuned on the carrier frequency o, = a ‘ 
and Jue 
20 
(38) 


where 2w is the total bandwidth between half power points. The current i(t) 
is the coaplex driving function, and its expression, if we assume K, = 1 and 
Ao RES e 1, follows from (6) 


jl w teeta ured] 


T(t) »« e+? (39) 


a rr SS PA CS SS PNP a PI ay TT ay LEST SO OPIN TP PETE WEA I Pt EET idee. | ee ee 
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From (2) we have 


P P 
i 1 = 
Ap) = 36 =" FE FF) 
a 

where 

P 

1 eee 1 j re i : Bese 1 

?, oRc = (.) 4 R22 .) Jie” 


Introducing now the Q@ of the circuit, we have 


i @, tr ~ = 24 2RC 
Py 1 1 
oc - w+ jo, Ap) = de (gee + gee) 
P, 1 2 
With reference to fomulae (25), (26) and (27), we recognize that 
B 
1 1 
t= lle sg oo. ys myer. 6, ~ =% 
Y § 
) a d a 0 ’ 1 = + © t- = ’ 1 = + © t 
1 2 v5 ° vy $0 


Bq. (17), taking into account (a4), (42), (43), becomes 


56 
= ee -wt+ SY 2 
Hy) « mite S790) «Be HOD f - i. 1 
sl 
we+—2e 57 3, 
+21, (A,) de » ome | 


(40) 


(41) 


(42) 


(43) 


(44) 
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where, recalling Eqs. (31) and (32), 
v Sat 4 
A ao (we 3%) = =m Wt = (ee, t+ 2) (45) 


We now want to investigate the magnitude of the term Containing 


rat Aa} in Eq. (44), and compare it with the mgnitude of the term containing 
TQ). in order to see if some simplication is possible. 
In the following we shall assume that 
20 
- (46) 
Let's start investigating the term containing 1,(A,). Ye have 
2 
Iagt = [ OG = wy? + (ant + ¥)? | 
2 ° 
02 2 2. (Me V2 bit 
=2) + (wt) + (20, t) +) > [ (47) 


Comparing (47) and (46) we can conclude that cliotar ais From Bq. (36) we get, 


for the magnitude of the term containing Tha o) 


a er 


fo} Vv ®, 
sar owt-J(20, t 


Let's now investigate the term, containing I Ay) We have, letting Z= §- 


Te EE Ee TS SEPT ETETE PPP PUTT EL She FOP | ER PLS PT Se * 
ny 
m4 
Tw 
oo 
Y 
nm 
~~ 
Fe. 
@ 
| 
+ 
e. 
=< 
Tw 


s w-j 2 
-wt + jv — -% 
oy go eo 
~wt -(§ -3 ¥) ~wt -( € -5 ¥) 
~2e% e ar 2 *, e a 
; rt ale lst el §-3% s | 


OPW ek eh Em TS te tes The? eerie Tee pee os ese 


A a a ee 
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Since in virtue of (36) 
aes 


we can conclude that, in Eq. (44), the term containing I,(A.) can be neglected 
in comparison with the term containing I,(A,). With thif sfmplitication, and 

keeping into account (38), we oasily obtain the following expressions for the 

real and imaginary part of v(t): 


1} [2 wt v_ vy 2 
Reo ES (Try cos Tea vin +d | cone, t 


x 


2" 4 x 
+S (lay eine +I. 5 cos re sinw ¢ (50) 


a 


= a Loe ¥ ¥ 
b ¢ ds rg ( re cosy Mo sin ©) cos wt 


+ i Y (“Toa sin z+ Tho Cos *) + z | sin w, % (51) 
or letting 

A = + f eo (Thy cos . - re sin re + , | (52) 

Boe 5 few (Tia vin + 1 cos 7a (53) 


2 = Acosw, t+ 8 sing, t (54) 

X = -Boosw, t+Asinw t (55) 
or 

Bw Va> +B" cos (o, t - $) Ys tan” 2 (56) 

x = Va? 43° cos (ow, t + 5 - $) (57) 


Except for the constant K o, we can therefore write 


v 


ad 
~waox at (58) 


a ES SS a A a NEY ee Te TE) SES | SSR oon i 
j 
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or 
BY - Ba! 
vq(t) ——— Be 5 awcrk (59) 
os tae 
a2 % 
Ay =Xwha2te Ld fr4¢9) cos S = ICA) vin ¥] 4 (60) 
ie a 
2° eewB=2%o » [F,in) sin + 1 (a) cos | (61) 


and taking into account the first cf Bq. (45) 


2 
20k! « -L a1 (62) 
ot leu 
2cB! 2 -B+2 — : (63) 
leu 
where we have 
us pt (64) 


Section V_ - Asymptotic Formulae 


It is now interesting to investigate the form that Bq. (59) takes 
when w—?00 or w—-+9%. The first case obviously corresponds to the case of 
@ pure resistance across the current generator (see Fig. 3). We expect there~ 
fore to find Ya proportional to the modulating function (5). The second case 
corresponds to a lossless parallel tunec circuit across the current generator. 
However, in this case the form of va is not so obvious as in the first case. 


To find the asymptotic formula for v, when w—>o we use the asymp- 
totic expression (36). In order to get a remainder for the quantities 4), B). 
A', B' that vanishes at least as ww, we have to take two (or more) terms of 
the expansion (36). Proceeding in this way we get 
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pret See yaaa ta 
oe war ie 2b bass + 0 ch) (66) 
2Ca! = +4 mae +0) 67) 
208 « + pe oe + 0 (4b) (68) 


introducing these expressions in Eq. (59), and keeping only the term of lowest 
order in p/w, we get 


lim v,(t) = 2h 2 5 (69) | 
w— 700 leu 


the right hand side of Eq. (69) is simply d/dt (2 tan wt) i.e. the modulating 
function 8(t). It shonld be noted that the asymptotic expression (69) could have 
aleo deen obtained by introducing the first term of Eqs. (36) in (60) ana (61) 
and then evaluating © A’ and C B’. Such a yrocedure is equivalent to differ- 
entiating the asymptotic expansion (36). In this case it is easy to see that 
one term only of the asymptotic expression leads to a remainder which vanishes 
at least as b/w. To find the asyuptotic formule for w/p—30, we use the ceries 
expansion (33). Proceeding in this way, and keeping only the lowest order in vw, 
we get 


lim v, = - owK—% = (t) (70) 


Section VI - Numerical Computations 


Using Eqs. (60), (61), (62), (65) and (59) numerical computations 

have been carried out for - 4< u = pt < 4 for tp follbwingwalues of w/p =-50, 

10, 4, 2, 1, OU, 0.1, 0.01, ——» 0. The results of these calculations have been 

recorded in Pigs. X8I-13106, MRI-13107, MRI-13108, MRI-13109, MRI-13110, MRI-13111, 

MEI-13112, MBI-13113, NRI-13114 and MEI-i3115, where v, 1s plotted sgminst u, 

and we have let vA - v4/%s of?) 1 oo 6) = ok as if4s apparent from Iq. (70). 
. * 


i 


re Te Ye 


a PY TE ET i a ee ae era eT 
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The naxinug positive value of vi bas deen normalised to 1. 1% get the actual 
value of v,, 1. @. the valus (59$ divided by 2y K, the ordinates of the diagram 
have to be multiplied by a constant EH whose value is indicated in each diagran. 


The numerical evaluation of the integral (3C) has been carried out 
by means of the expansions (35) or (36). When possible the values tabulated by 
the Computation Labcrstory, 150 Nassau Straet, New York, N.Y. have been used. 
However, the available Tables are, at present, rather limited, and in particular 
the values corrssponding to negative real part of the independent variable, have 
not been tabulated. 


The resuits of the computation for I.(A) end I ALS, have been recorded 
in Figs. MRI-13116, MRI-13117, MRI-13118, MRI-13119, mptt131120 and MRI-13121. 
In these diagrams the points marked with an arrow refer ic the values tabulated 
by the National Laboratory. Ths other values have been computed by us. 

In the computation we have let 


AaexeJsy 


where 
w w 
@- -_ u - = 
x =z rn ¥ 2 i 
Conclusions 


From the spt of Fige. MRI-13106 to MRI-13115 it 1s apparent that the 
tine dependence of v, 1s very similar to that of S(t) if w/p >1. This result 
ig interesting since it shows that, even with a relatively narrow passband, the 
reproduction of the priginal signal is still good. To appreciate this fact con- 
sider a pulse of 10°~° sec. and assume for duration of the pulge the definition 
given in Section IZ. From Eq. (21) therefore we get » = 6.10%, so that in order 
to get a good response it is necessary to assuze >6.10°, or the total band- 
width in cycles greater than approximately 210%o/seo. Figures MRI-13111 and 
MBI-13112 show, however, that it is not safe to operate in the neighborhood of 
w/p = 1, since for instance for w/p = 0.4 the response is very distorted. 


It has to be emphasized that the numerical results obtained refer 
only to the case of a =H. m/s = 1. It would be interesting to carry on calcu- 
lations for other values of q and eventually ootain formulae for non-integral 
walues (<1) of q. However, it seems that the formulae become very cumbersome 
even for values of d= 2 or 3. It is also doubtful that it would be feasible 
to obtain formulae for q 73 independent of @, + 
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